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Aims. We probe the dependence of the electron-to-proton mass ratio, /i = m c /m p , on the ambient matter density by 
means of radio astronomical observations. 

Methods. The ammonia method, which has been proposed to explore the electron-to-proton mass ratio, is applied 
to nearby dark clouds in the Milky Way. This ratio, which is measured in different physical environments of high 

O (terrestrial) and low (interstellar) densities of baryonic matter is supposed to vary in chameleon-like scalar field models, 
which predict strong dependences of both masses and coupling constant on the local matter density. High resolution 
spectral observations of molecular cores in lines of NH3 (J, K) — (1, 1), HC3N J = 2 — 1, and N2H + J = 1 — were 
performed at three radio telescopes to measure the radial velocity offsets, AV = V ro t — Vinv, between the inversion 
transition of NH3 (1,1) and the rotational transitions of other molecules with different sensitivities to the parameter 

A/i//X = (/i obs — ^lab)/^lab- 

Results. The measured values of AV exhibit a statistically significant velocity offset of 23 ± 4 sta t ± 3 sys m s _1 . When 
interpreted in terms of the electron-to-proton mass ratio variation, this infers that A/1//1 = (2.2±0.4 sta t ±0.3 sys ) x 10 -8 . 
' If only a conservative upper bound is considered, then the maximum offset between ammonia and the other molecules 

I is |AV| < 30 m s _1 . This provides the most accurate reference point at z = for Afi/fi of |A/i//i| < 3 x 10 -8 . 
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£^ . 1. Introduction cause of its coupling to matter that strongly depends on 
J^J the ambient matter density. In this way, this scalar field 
; Testing the variability of dimensionless physical constants may evade local tests of the equivalence principle and fifth 
■ I . is an important topic in contemporary laboratory and as- force experiments since the range of the scalar-mediated 
trophysical experiments. The physical constants are not fifth force for the terrestrial matter densities is too small 
_ supposed to vary in the Standard Model (SM) of parti- to be detected. This is not the case for space-based tests, 
CTS . cle physics but do vary quite naturally in grand unification where the matter density is considerably lower, an effec- 
■ theories, multidimensional theories, and whenever there is tive mass of the sca i ar fi e ld is negligible, and an effective 
a coupling between light scalar fields and baryonic matter. range for the sca l ar _ m ediated force is large. The present pa- 
in particular, light scalar fields have been widely discussed per dea ls with one such astronomical test based on spectral 
in the context of dark energy, since they provide negative observations of molecular clouds in the Milky Way disk, 
pressure that may be responsible for the cosmic accelera- Additional tests employing polarization of the light from 
tion detected at low redshifts, z < 1 (Caldwell et al. 1998; the stars and a modification of the Sunyaev-Zel'dovich ef- 
Peebles & Rata 2003). If this scalar field does exist, then a fect in the cosmic microwave background due to a coupling 
question arises: why has it not been detected in local tests between a chameleon-like scalar field and photons are de- 
of the equivalence principle or fifth force searches? A so- scribe d, respectively, in Burrage et al. (2009) and Davis et 
lution was suggested using the so-called 'chameleon' mod- a j (2009) 
els (Khoury & Weltman 2004; Brax et al. 2004; Mota & 

Shaw 2007). These models assume that a light scalar field Astronomical spectroscopy can probe the physical con- 
acquires both an effective potential and effective mass be- stants that describe atomic and molecular discrete spectra: 

the fine-structure constant, a = e 2 /(frc), and the electron- 
* Based on observations obtained with the Medicina 32-m tele- ^-proton mass ratio, /i = m c /m p . In GUTs, these con- 
scope operated by INAF - Istituto di Radioastronomia, the 100- stants mediate the strength of fundamental forces: a is the 
m telescope of the Max-Planck Institut fur Radioastronomie at coupling constant of the electromagnetic interaction, m e is 
Effelsberg, and the Nobeyama Radio Observatory 45-m tele- related to the vacuum expectation value of the Higgs field, 
scope of the National Astronomical Observatory of Japan. namely to the scale of the weak nuclear force, and m p is 
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proportional to the quantum chromodynamics scale. We 
note that the predicted variabilities of the fine-structure 
constant and the electron-to-proton mass ratio are not in- 
dependent and that the variations in p may exceed those 
of a (Calmet & Fritzsch 2002; Langacker et al. 2002; Dine 
et al. 2003; Flambaum et al. 2004). Thus, a hypothetical 
variation in p is expected to be easier to detect than that 
in a. 

Despite many efforts, the variability in a and p has 
never been detected. Data obtained from high precision 
frequency measurements in laboratory experiments with 
atomic clocks and from astronomical observations pro- 
vide only upper limits. For example, laboratory experi- 
ments have delivered the following results: fractional tem- 
poral variations in p are restricted to a level of p/p = 
(3.8 ± 5.6) x 10~ 14 yr- 1 (Shclkovnikov et al. 2008), and 
p./ 'p = (1.6 ±1.7) x 10~ 15 yr" 1 (Blatt et al. 2008), whereas 
the current level for a is a/a — (—1.6 ± 2.3) x 10~ 17 yr _1 
(Rosenband et al. 2008). Here Ap and Aa are the rela- 
tive changes between the values measured at two different 
epochs. 

The most stringent upper limits to Aa/a and Ap/ p 
obtained from astronomical observations of extragalactic 
objects are restricted by a few units of ppm (lppm = 10~ 6 ). 
Here Aa/a = (a' — a)/ a, where a, a' denote the values of 
the fine-structure constant in the laboratory and the spe- 
cific absorption/emission line system of a Galactic or extra- 
galactic object (the same definition is applied to Ap/p). 
However, there are claims of a variability in both a and p 
at the 5cr and 4cr confidence levels, respectively, although 
they are in contrast to other null results and the whole issue 
is highly controversial. The current observational status is 
the following. 

From the measurements of the relative radial velocity 
shifts of different absorption lines (e.g., Mgn, Sin, Fen, 
Znn) of 143 QSO absorption systems obtained with the 
Keck/HIRES in the redshift range 0.2 < z < 4.2, Murphy 
et al. (2004) claimed that Aa/a = -5.7 ± 1.1 ppm. On the 
other hand, Chand et al. (2004) using VLT/UVES spectra 
of bright QSOs analyzed 23 absorption systems, which are 
not in common with those studied by the Keck telescope, 
and failed to reproduce Murphy et al.'s result. At first, 
they claimed a very stringent limit of Aa/a — — 0.6 ± 0.6 
ppm, but Murphy et al. (2008a) demonstrated convincingly 
that the error estimate of Chand et al. (2004) is underesti- 
mated. When these errors are properly accounted for, the 
new weighted mean becomes Aa/a = — 6.4 ± 3.6 ppm. The 
new value has a scatter that is larger than the quoted er- 
rors implying that there are systematic errors that are com- 
parable or even larger than the statistical ones. However, 
Srianand et al. (2008) showed that excluding two systems 
that deviate from the mean at the 3a level when reanalyzing 
VLT/UVES systems leads to Aa/a = 0.1 ± 1.5 ppm. We 
note that the only claimed non-zero measurement of Aa/a 
corresponds to an uncertainty in the wavelength scale cal- 
ibration smaller 100 m s" 1 . This precise calibration of the 
Keck/HIRES spectra is difficult to achieve as shown by 
Griest et al. (2010). Using iodine exposures to calibrate the 
normal Th-Ar Keck data pipeline output, they found ab- 
solute wavelength offsets of 500 m s _1 to 1000 m s" 1 with 
drifts of more than 500 m s _1 over a single night, and drifts 
of nearly 2000 m s _1 over several nights. Their conclusion is 
that these systematic uncertainties make it difficult to use 
Keck/HIRES QSO spectra to constrain the change in the 



fine structure constant at the 10 -6 level. The most strin- 
gent constraint on |Aa/a| < 2 ppm was found from the 
Fen lines at z = 1.15 towards the bright quasar HE 0515- 
4414 (Quast et al. 2004; Levshakov et al. 2005; Molaro et 
al. 2008). 

Bounds on p variations are most effectively obtained 
from observations of the Werner and Lyman series of the 
molecular hydrogen H2 in damped Ly-a systems (DLA). 
The electron- vibro-rotational transitions have different de- 
pendences on the reduced mass and can be used to con- 
strain the variability of p (Thompson 1975; Varshalovich 
& Levshakov 1993). The measurements of p rely on the 
same H2 systems observed with VLT/UVES and claim 
a variability of Ap/ p = 24 ± 6 ppm (Reinhold et al. 
2006) or no variability with —15 ppm < Ap/p < 57 
ppm (Levshakov et al. 2002), \Ap/p\ < 49 ppm (Wendt & 
Reimers 2008), Ap/p = 2.6 ± 3.0 ppm (King et al. 2008), 
and Ap/p = — 7± 8 ppm (Thompson et al. 2009). Radio 
astronomical observations place constraints on p- variations 
of \Ap/p\ < 1.8 ppm at redshift z = 0.68 (Murphy et al. 
2008b), and \Ap/p\ < 0.6 ppm at z = 0.89 (Henkel et al. 
2009). 

One should keep in mind, however, that laboratory ex- 
periments and spectra of extragalactic objects probe very 
different time-scales and different regions of the universe, 
and the connection between them is quite model dependent 
(Mota & Barrow 2004). 

Astronomical measurements of the dimensionless con- 
stants are based on the comparison of the line centers in the 
absorption/emission spectra of astronomical objects and 
the corresponding laboratory values. To distinguish the line 
shifts due to the radial motion of the object from those 
caused by the variability of constants, lines with different 
sensitivity coefficients, Q, to the variations in p and/or a 
should be employed. Unfortunately, the observable opti- 
cal lines have very close sensitivity coefficients with differ- 
ences of |AQ| that do not exceed 0.05. Combined with the 
resolving power of spectrographs at modern optical tele- 
scopes, this small difference leads to an upper limit on 
Aa/a and Ap/p of ~1 ppm, the optimal value achiev- 
able in observations of extragalactic objects with existing 
facilities. However, one can probe the values of fundamental 
constants at a considerably more accurate level if, firstly, 
nearby objects in the Milky Way are observed and, sec- 
ondly, spectral lines from other frequency ranges, infrared 
and radio, are employed. This statement is based on the 
following considerations. 

As mentioned above, the scalar field models suggest 
a coupling between the scalar fields and baryonic matter. 
This coupling results in a functional dependence of p and a 
on p, the local matter density (Olive & Pospelov 2008). This 
coupling alone may prevent a positive detection of the vari- 
ations in p and other dimensionless constants in laboratory 
studies, since they are performed under the same terres- 
trial conditions and the effective range of the scalar-ficld- 
mediated force is smaller than 1 mm at terrestrial matter 
densities (Olive & Pospelov 2008). On the other hand, the 
density in cold molecular clouds is only 10 3 — 10 5 particles 
per cm 3 , i.e., ~ 10 19 times lower than in terrestrial environ- 
ments. Because of this extremely large difference between 
the matter densities non-zero values of Ap/p and Aa/a 
are predicted. It is very important to compare the above 
ratio of ~ 10 19 with the differences in p between molecu- 
lar clouds themselves, to find that they are negligible, i.e., 
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at most one or two orders of magnitude for a given tracer. 
In the context of Afi/fi (or Aa/a) this means that mea- 
surements of all dense molecular clouds are identical irre- 
spective of their location in space and time (redshift). Thus, 
nearby Galactic molecular clouds can be observed to ensure 
a strong signal, i.e., that the line profiles can be centered 
with high accuracy. 

Molecular lines originating in these clouds are mainly 
observed in cm and mm radio bands. An advantage of ra- 
dio observations is that very narrow spectral lines (of line 
widths -100 m s" 1 ) arising in cold molecular cores can be 
observed with extremely high spectral resolution, FWHM 
~ 25 m s _1 . In optical observations of extragalactic ob- 
jects, we have in general broader sources, b > 1 km s _1 , and 
use lower resolution spectrographs, FWHM > 4 km s _1 
(Levshakov et al. 2007). Taking into account that the un- 
certainty in the line position is roughly l/10th of the pixel 
size, this infers a gain of about two orders of magnitude in 
the accuracy of the measurements of Aa/a and/or A/i/ /i 
if radio spectra of the local interstellar objects are used. 

Complementary to this, differences in the sensitivity co- 
efficients for lines from the microwave, far IR, and radio 
ranges are much larger than those from optical and UV 
ranges. For example, the difference between the sensitiv- 
ity coefficient of the (J,K) = (1,1) inversion transition of 
ammonia (NH3) and that of any rotational transition in 
another molecule is AQ = 4.46 — 1 = 3.46 (Flambaum & 
Kozlov 2007). When compared with UV transitions, this 
infers a gain of about 70 times in the sensitivity of the line 
positions to the change in /i. 

Our preliminary study of cold molecular clouds 
(Levshakov et al. 2008a, hereafter LMK) was based on pub- 
lished results of high quality spectral radio observations ob- 
tained with the 100-m Green Bank telescope (Rosolowsky 
et al. 2008; Rathborne et al. 2008) and the 45-m Nobeyama 
radio telescope (Sakai et al. 2008). Comparison of the rel- 
ative radial velocities of the NH3 (J,K) = (1, 1) inversion 
transition and CCS Jjv = 2i — lo rotational transition mea- 
sured a systematic shift of AV = V/ sr .(CCS) — VJ sr (NH 3 ) — 
60 m s _1 , which when interpreted in terms of the electron- 
to-proton mass ratio variation infers that A/1//1 — 6 x 10 -8 . 
We also noted that a similar offset between N2H 4 " and NH3 
measured by Pagani et al. (2009) in the cold dark cloud 
L183 could have a similar physical origin (Molaro et al. 
2009). 

In this paper, we present results of our own spec- 
tral observations of cold and compact molecular cores in 
the Taurus giant molecular complex obtained with the 
Medicina 32-m telescope, the Effelsberg 100-m telescope, 
and the Nobeyama 45-m telescope. 



2. The ammonia method 

Narrow molecular lines observed in cold dark clouds provide 
a sensitive spectroscopic tool to study relative shifts of the 
order of a few 10 ms" 1 between radial velocities of different 
molecular transitions. Among numerous molecules detected 
in the interstellar medium, ammonia (NH3) is of particu- 
lar interest because of the high sensitivity of the inversion 
frequency to a change in /1. The inversion vibrational mode 
of NH3 is described by a double-well potential, the first 
two vibrational levels lying below the barrier. The quantum 
mechanical tunneling splits these two levels into inversion 



doublets providing a transition frequency that falls in the 
microwave range (Ho & Townes 1983). 

For the ammonia isotopologue 15 ND3, van Veldhoven et 
al. (2004) first showed that the inversion frequency of the 
(J, K) — (1, 1) level varies a£] 



^) =5.6^ 
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i.e., the inversion transition is an order of magnitude more 
sensitive to /x-variation than molecular vibrational transi- 
tions, which scale as i5 viD — /i 1 / 2 . 

The sensitivity coefficient of the inversion transition 
(J,K) = (1,1) in NH 3 (v = 23.7 GHz) was calculated by 
Flambaum & Kozlov (2007) from the numerical integration 
of the Schrodinger equation for different values of /z, and 
from the analytical Wentzel-Kramers-Brillouin (WKB) ap- 
proximation of the inversion frequency. Both methods pro- 
vide similar results giving 



Av 

v 
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Here v and v are the frequencies corresponding to the lab- 
oratory value of /i and to an altered /1 in a low-density 
environment, respectively. 

The rotational frequency scales as E rot — /i and, thus, 



Av 
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By comparing the observed inversion frequency of NH3 
(1,1) with a suitable rotational frequency of another 
molecule produced co-spatially with ammonia, a limit on 
the spatial variation of [i can be determined. 

In radio astronomical observations, any frequency shift 
Av is related to the radial velocity shift AV r (V r is the 
line-of-sight projection of the velocity vector) 



AV r 



V r -Vo 



flab — fobs 
^lab 



(4) 



where Vq is the reference radial velocity, and fi a b, fobs 
are the laboratory and observed frequencies, respectively. 
Therefore, by comparing the apparent radial velocity Vj nv 
for the NH3 inversion transition with the apparent radial 
velocity V ro t for rotational lines originating in the same 
molecular cloud and moving with a radial velocity Vq with 
respect to the local standard of rest, we can find from 
Eqs. P-iD 



Ay 
A* 



0.289- 



V r 



0.289- 



AV 



(5) 



From Eq.([S]), we can estimate the limiting accuracy of 
Afi/ fi available from modern radio astronomical observa- 
tions. At high spectral resolution in the microwave range 
(FWHM — 30 m s _1 ), the errors in the molecular line 
position measurements are mainly restricted by the uncer- 
tainties in laboratory frequencies, e v — 0.1 — 1 kHz, which 
correspond to the Vi sr uncertainties of e v — 1 — 10 m s _1 . 
Taking into account that A/.t//i — 0.3AV/c, the expected 
accuracy of A////1 in a single measurement is about 10 -8 . 



1 The sign of the sensitivity coefficient Q = 5.6 was misprinted 
in van Veldhoven et al. (2004) , as noted by Flambaum & Kozlov 
(2007). 
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This level of accuracy can be achieved, however, only un- 
der ideal conditions since it assumes that molecules are 
identically distributed within the cloud and are observed 
simultaneously with the same receiver, beam size, system 
temperature, and velocity resolution. Violation of any of 
these conditions leads to random shifts of the line centers, 
which are referred to as the Doppler noise. The input of the 
Doppler noise to a putative A\ij \i signal can be reduced to 
some extent if the velocity shifts due to the inhomogeneous 
distribution of molecules and to instrumental imperfections 
are minimized. For our purposes, we express the velocity 
offset AV in Eq.© as the sum of two components 

AV = AV^ + AV n , (6) 

where AV^ is the shift due to //-variation and AV n is the 
Doppler noise. We assume that the noise component has 
zero mean (AV n ) = km s _1 and a finite variance. The 
signal AV M can then be estimated statistically by averaging 
over a data sample 

(AV) = (AVfj,), Var(AV) = Var(AV fi ) + Var(AV n ). (7) 

From this equation it is seen that the measurability of the 
signal (AV/j,) depends critically on the value of the Doppler 
noise. The Doppler noise Var(AV n ) can be reduced by an 
appropriate choice of molecular lines and by special criteria 
applied to the selection of targets. 

2.1. Molecules appropriate for Afi/ fi measurements 

The observed molecular transitions should as far as possible 
share the same volume elements to have similar Doppler ve- 
locity shifts. The ammonia inversion transitions are usually 
detected in dense molecular cores (n > 10 4 cm -3 ), which 
are represented in the Milky Way disk by a large variety 
of types and physical properties (Di Francesco et al. 2007). 
Mapping of the dense molecular cores in different molecu- 
lar lines shows that there is a good correlation between the 
ammonia NH 3 , N 2 H+, and HC 3 N distributions (e.g., Fuller 
& Myers 1993; Hotzel et al. 2004; Tafalla et al. 2004; Pagani 
et al. 2009). However, in some clouds NH 3 is not traced by 
HC3N. The most striking case is the dark cloud TMC-1, 
where peaks of line emission are offset by 7 arcmin (Olano 
et al. 1988). 

In general, N-bearing molecules such as NH3 and N2H 4 " 
trace the inner cores, where the density approaches 10 5 
cm~ 3 . At the same time, the carbon-chain molecules dis- 
appear from the gas-phase because of freeze-out onto dust 
grains (Tafalla et al. 2004). HC 3 N, as well as other C- 
bearing molecules, are usually distributed in the outer parts 
of the cores. The mutual distribution of NH 3 and HC 3 N is 
affected by chemical differentiation in the process of the 
dynamical evolution of the core. HC 3 N is abundant in the 
early evolutional stage of the star-forming regions (Lee et 
al. 1996), when the fractional abundance of HC 3 N remains 
almost constant and the spatial distributions of the N- and 
C-bearing molecules match each other quite well (Suzuki 
et al. 1992). As the gas density increases and at dust tem- 
peratures Td < 20 K, these distributions diverge because of 
adsorption of the heavy molecules (e.g., HC 3 N) from the gas 
phase onto grain mantles (Flower et al. 2006). In the later 
stages of the protostellar collapse, carbon chain molecules 
are destroyed by high velocity outflows and radiation from 
protostars, whereas the same processes favor the desorption 
of ammonia from dust grains (Suzuki et al. 1992). 



The chemical differentiation and velocity gradients 
within the molecular core are the main sources of the un- 
avoidable Doppler noise in Eq.©. Additional scatter in 
the AV n values is caused by the different optical depths of 
the hyperfine structure transitions. However, all of these ef- 
fects lead to the radial velocity shifts detected between NH 3 
and other molecules (HC 3 N, N2H + ), which may be random 
from core to core. Thus, being averaged over a large sam- 
ple of targets, the Doppler noise component AV n in Eq. ^ 
should be canceled out. 



3. Observations 

For relative velocity measurements, the molecular lines of 
NH 3 , HC 3 N, and N2H" 1 " were chosen. From the published 
data, we selected 41 molecular cores with compact mor- 
phology, i.e., the cores whose geometrical structure can be 
represented by a central region of nearly constant density 
and a surrounding envelope with a density that decreases 
as a power law. For the dark clouds in the Taurus molecular 
complex, a typical size of the central core is ro ~ 20" — 60" 
(e.g., Daniel et al. 2007). The selection was based on nar- 
row and sufficiently strong molecular lines that correspond 
to individual hyperfine transitions: NH 3 (J,K) = (1,1), 
HC 3 N J = 2 - 1, and N 2 H+ J = 1 — 0. This suggests that 
the selected cores have low kinetic temperatures (Tki n — 10 
K) and are located at distances below 250 pc. The kinetic 
temperature is known to be surprisingly uniform in dark 
clouds (e.g., Dickman 1975; Walmsley & Ungerechts 1983). 
The list of sources is given in Table [T] Observations were 
performed with radio telescopes in Medicina, Effelsberg, 
and Nobeyama between November 2008 and April 2009. 

Medicina. Observations at the Medicina 32-m tele- 
scope0 were carried out on November 24-28, 2008. 
Both available digital spectrometers ARCOS (ARcetri 
Correlation Spectrometer) and MSpecO (high resolution 
digital spectrometer) were used with channel separations 
of 4.88 kHz and 2 kHz, respectively. For ARCOS, this cor- 
responds to 62 m s _1 at the position of the ammonia in- 
version transition (23 GHz) and 80 m s _1 at the rotational 
HC 3 N (2-1) line (18 GHz). For MSpecO, it is 25 m s" 1 
and 32 m s _1 at the corresponding frequencies. Only a few 
brightest objects (marked by a symbol m in Table [T]) were 
observed. The Medicina 32-m telescope angular resolution 
is -1.6' at 23 GHz and -2.1' at 18 GHz. The pointing 
accuracy was superior to 25". Spectra were taken in a po- 
sition switching mode with a typical integration time of 5 
min for both ON- and OFF-source scans. The OFF posi- 
tion was taken to be approximately five beam widths to 
the west of the source position. Typically, 10-20 ON/OFF 
pairs were taken, depending on the frequency and source 
flux. Unfortunately, because of poor weather conditions not 
all observations studied both NH 3 and HC 3 N transitions. 
Standard data reduction was performed using the CLASS 
reduction packag^E 

Effelsberg. The (J ,K) = (1,1) inversion line of ammo- 
nia (NH 3 ) and the J=2-l rotation line of cyanoacetylene 
(HC 3 N) were also observed with the 100-m telescope at 



2 The 32-m VLBI antenna at Medicina is operated by the 
INAF-Istituto di Radioastronomia in Bologna. 

3 http: / /www. iram.fr/IRAMFR/ GILDAS / 
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Effelsbergjj] on February 20-22, 2009. The corresponding 
targets are marked by symbol e in Table [TJ The lines were 
measured with a K-band HEMT (high electron mobility 
transistor) dual channel receiver, yielding spectra with an- 
gular resolutions of 40" (NH 3 ) and 50" (HC 3 N) in two 
orthogonally oriented linear polarizations. Averaging the 
emission from both channels, typical system temperatures 
are 100-150 K for NH 3 and 80-100 K for HC 3 N on a main 
beam brightness temperature scale. 

The measurements were carried out in frequency switch- 
ing mode using a frequency throw of ^5 MHz. The backend 
was an FFTS (fast fourier transform spectrometer), oper- 
ated with its minimum bandwidth of 20 MHz, providing 
simultaneously 16 384 channels for each polarization. The 
resulting channel widths are 15.4 and 20.1 m s _1 for NH 3 
and HC 3 N, respectively. We note, however, that the true 
velocity resolution is about twice as large. 

Observations started by measuring the continuum emis- 
sion of calibration sources (NGC 7027, W3(OH), 3C286) 
and continued by performing pointing measurements to- 
ward a source close to the spectroscopic target. Spectral 
line measurements were interspersed with pointing mea- 
surements at least once per hour. The calibration is esti- 
mated to be accurate to ±15% and the pointing accuracy 
to be superior to 10 arcsec. The CLASS reduction package 
was used for standard data reduction. 

Nobeyama. The NH 3 (J, K) = (1,1) line at 23 GHz 
and the N2H + J =1-0 line at 93 GHz were observed with 
the Nobeyama Radio Observatory (NRO) 45-m telescop^ 
on April 8-10, 2009. We used a low-noise HEMT receiver, 
H22, for the NH 3 observations and the two sideband- 
separating SIS (Superconductor-Insulator-Superconductor) 
receiver, TlOO (Nakajima et al. 2008), for the N 2 H+ obser- 
vations. Both of them are dual polarization receivers. We 
observed two polarizations simultaneously. The 45-m radio 
telescope angular resolution is about 73 and 17 arcsec at 
23 and 93 GHz, respectively, and the main beam efficiency 
is 0.84 at 23 GHz, and 0.53 at 93 GHz. Autocorrelators 
were employed as a backend with bandwidth and channel 
separation of 4 MHz and 4.375 kHz, respectively. This cor- 
responds to channel widths of 57 ms" 1 at 23 GHz, and 14 
m s- 1 at 93 GHz. 

The telescope pointing was checked by observing nearby 
SiO maser sources every 1-2 hr. The pointing accuracy was 
< 5". The line intensities were calibrated using the chopper 
wheel method, and the observations were carried out in po- 
sition switching mode. The data reduction was performed 
partly with the AlPS-based software package NewStar de- 
veloped at NRO. 

4. Data analysis 

4.1. Vis,, calculation 

The radial velocities, VJ sr , are determined from the spec- 
tral line analysis. Each individual exposure was first visu- 
ally analyzed and corrupted data were excluded. Individual 
exposures were then coadded to increase the signal-to-noisc 



4 The 100-m telescope at Effelsberg/Germany is operated by 
the Max-Planck-Institut fur Radioastronomie on behalf of the 
Max-Planck-Gesellschaft (MPG). 

5 The 45-m radio telescope is operated by Nobeyama Radio 
Observatory, a brach of the National Astronomical Observatory 
of Japan. 



(S/N) ratio. Before coadding, a baseline was removed from 
each spectrum. To define the baseline, spectral intervals 
without emission lines and/or noise spikes were selected 
and the mean signal Tj along with its rms uncertainty cii 
were calculated for each interval. A set of pairs {Tj, <ii\ was 
used to find a baseline (regression line) by minimizing \ 2 - 
This baseline was typically linear but was in some cases 
quadratic or cubic. Since individual rms uncertainties <Ji 
were of the same order of magnitude, their mean value a 
was assigned to the whole spectrum. The resulting spectra 
were coadded with weights that are inversionally propor- 
tional to their variances, a 2 . 

The line parameters such as the total optical depth in 
the transition, t, the radial velocity, V/ sr , the line broad- 
ening Doppler parameter, 6, and the amplitude, A, were 
estimated by fitting a one-component Gaussian model to 
the observed spectra. The model was defined by 



T(v) 



where 



t(v) 



A' 



-*(«) 



k 

E 



a, exp 



(v 



V lsr ) 2 



IP 



(8) 



(9) 



where a^, Vi is, respectively, the relative intensity of the ith 
hypcrfinc component and its velocity separation from the 
reference frequency. The sum in Eq. ([9]) runs over the k 
hyperfine structure (hfs) components of the transition. The 
physical parameters <2j, and u, for the NH 3 (1,1), HC 3 N 
(2-1), and N 2 H+ (1-0) transitions are listed in TablesEJH 
and|H respectively. 

The fitting parameters {A, r, b, Vi sr } were determined 
by means of a ^-minimization procedure. For optically thin 
transitions, Eq.® transforms into 



T(v) = A-t(v), 



(10) 



which prevents the independent estimation of A and r. In 
this case, the model parameters were {A ■ t, b, Vi sr }. 

Since we are mostly interested in the model parameters 
Vi sr and b, their values are listed in Tables [S] and [51 The la 
errors of Vi sr and b were estimated from the diagonal ele- 
ments of the covariance matrix calculated for the minimum 
of x 2 - The error in Vi sr was also estimated independently 
by the A% 2 method (e.g., Press et al. 1992) to control both 
results. When two estimates differed, the larger error was 
adopted. 

4.2. Uncertainties in the rest-frequencies 

Molecular lines observed in cold dark clouds of T^ n ~ 10K 
can be extremely narrow with a thermal broadening of 
Vth 100 m s _1 for NH 3 or even lower for heavier 
molecules. They provide a sensitive spectroscopic mea- 
surement of shallow velocity gradients in molecular clouds 
(Lapinov 2006). At high spectral resolution (FWHM ~ 30 
m s" 1 ), the uncertainty in the line position measurement 
from radio astronomical observations can be as small as 
1 m s" 1 which is comparable to the precision available 
for laboratory rest-frame frequencies of the NH 3 (J, A") = 
(1, 1) transition (Kukolich 1967; Hougen 1972). The physi- 
cal parameters of the observed hfs components of NH 3 are 
listed in Table [5J Column 1 indicates the group numbers 
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shown in the upper panel of Fig. [5] The errors in the line 
positions given in parentheses in Col. 7 are smaller than 1 
m s _1 with a mean of e v = 0.58 m s _1 . However, labora- 
tory uncertainties for other molecules used in the present 
observations are significantly larger. 

The analysis of all available laboratory data about 
HC3N carried out by Miiller et al. (2005) and independently 
by Lapinov (2008, private comm.) shows good agreement, 
the latter results being of slightly higher precision {e v ~ 2.8 
m s _1 ). The physical parameters of the hfs components of 
the J = 2 — 1 transition are presented in Table [3J The 
line identification numbers from Col. 1 are depicted in the 
middle panel of Fig. [2] 

The rest-frame frequencies for the third molecule N2H 4 " 
(Table @| were taken from the Cologne Database for 
Molecular Spectroscopy (CDMS) described in Miiller et 
al. (2005). The group numbers from Table HI Col. 1 are 
also shown in Fig. [51 lower panel. The CDMS data have a 
factor of two smaller errors than the hyperfine frequencies 
estimated from observations of the molecular core L1512 
in both N 2 H+ (93 GHz) and C 3 H 2 (85 GHz). It is again 
assumed that these molecules are co-spatially distributed 
(Caselli et al. 1995). The error in e v ~13.5 m s" 1 (~ 4 
kHz) reported in Table 0] corresponds to the CDMS data. 
The mean offset between the data from Caselli et al. and 
the CDMS is AV = V Case m - Vcdms = 12 m s" 1 . 

We do not correct the N2H" 1 " frequencies for a 40 m s _1 
offset as by Pagani et al. (2009), who used the rest- frame 
frequencies of Caselli et al.. The offset of 40 m s^ 1 , or 28 
m s _1 when adopting the CDMS rest frequencies (Molaro 
et al. 2009), was determined by comparing the inversion 
NH3 (1,1) and rotational N 2 H + (1 — 0) transitions observed 
in the target L183; this correction is justifable only when 
Afj,/fj, = 0. Reliable calibrations of the hyperfine transition 
frequencies of N2H 4 " can be performed only by a high pre- 
cision laboratory measurement. 

We note in passing that the estimate of the mean radial 
velocity based on the simultaneous fitting of n hyperfine 
transitions should be more precise than that based on a sin- 
gle line. However, the improvement is not as high as 
because the relative positions of the individual hfs transi- 
tions are correlated. If we consider n velocity differences 
{v\ ~vq, V2 — VQ, . . ., v n — Vq}, where vq is the reference ve- 
locity, then it is easy to show that the correlation coefficient 
Kij between two of them (i ^ j) is given by 



In cases of equal accuracy, the weight u>i = 1/n for each i. 
Then e sys is approximately equal to 



1 



(1 + *?)(! + *?) 



(11) 



where Si — a Vi / a Vo and s$ = a Vj / a Vo . Taking into account 
that the laboratory errors in the hyperfine line positions, 
a Vi , are almost equal, we have Kij — n rs 1/2. 

The covariance matrix Cov{vi — Vo,Vj — vq) contains n 
diagonal terms a 2 and n(n — 1) non-diagonal terms kct 2 , 
where a 2 represents the variance in a single measurement. 
The error in the mean radial velocity caused by the labo- 
ratory uncertainties (referred to as £ sys hereafter) can be 
calculated as described by, e.g., Stuart & Ord (1994) 



OJiUJjCov(vi - V , Vj - V ) 

i=1 j=l 



1/2 



— y/n + n(n — 1)k 



a\/K 



(13) 



(12) 



Thus, the gain factor, s/k, is only about 0.7 for the 
molecules in question. 

5. Results 

We first consider the entire sample of sources listed in 
Table [T] By applying the same computational procedure to 
all available spectra from the three radio telescopes, we find 
that not all of the molecular profiles can be described ade- 
quately with a single-component Gaussian model. In total, 
we measured n — 55 molecular pairs. The corresponding ve- 
locity offsets AV are shown in Fig.[T0]by three types of sym- 
bols: filled squares and circles (NH3/HC3N pairs), and open 
circles (NH3/N2H+ pairs) are, respectively, data points 
from the 32-m, 100-m, and 45-m telescopes. Some molecu- 
lar cores (L183, L1495, TMC-1C) were partially mapped. 
The offset coordinates in arcsec with respect to the source 
positions from Table [T] are shown in parentheses. Two ob- 
jects (L183 and TMC-1C) were observed at the 32-m tele- 
scope with different spectrometers ARCOS and MSpecO, 
as indicated in the figure. Five sources, TMC-1/HC3N 
(Medicina), TMC-2, L1521F, and L1544 (Effelsberg) , and 
TMC-2 (Nobeyama) from Table [T] exhibit asymmetric pro- 
files that are not consistent with a simple Gaussian model 
(we do not show them in Fig. [TU)) . A considerable frac- 
tion of the molecular pairs from the total sample exhibit 
non-thermal motion, i.e., the Doppler broadening parame- 
ter 6(NH 3 ) < 6(HC 3 N, N 2 H+). 

Several sources were observed at different radio tele- 
scopes. Since these data points have different systematic 
errors, we treat them as 'independent' measurements in the 
following statistical estimates. 

The weighted mean (weights inversionally proportional 
to the variances) of the ensemble of n = 55 AV values 
is (AV) W = 27.4 ± 4.4 m s _1 , the scale (standard devia- 
tion) is 32.6 m s _1 , and the median is 17 m s _1 . We also 
used a robust redescending M-estimate for the mean and 
the normalized median absolute deviation (1.483-MAD) for 
the scale. These statistics work well for inhomogeneous data 
sets with outliers and deviations from normality (corre- 
sponding formulae are given in Appendix). The M-estimate 
infers that (AV% = 14.1±4.0 m s" 1 (scale 29.6 m s" 1 ). A 
poor concordance between three mean estimates ((AV) W , 
(AV) M , and the median) is caused by large systematic shifts 
and, as a result, by 'heavy tails' of the probability distri- 
bution function. The scatter in the points reflects effects 
related to the gas kinematics and the chemical segregation 
of one molecule with respect to the other. 

An additional decrease in the noise component 
Var(AV n ) in Eq. (JJJ is possible if we select from the sam- 
ple of the observed targets the systems with 'simple' ge- 
ometry and internal kinematics. The ideal target would be 
a homogeneous spherical cloud where different molecules 
are co-spatially distributed (no chemical segregation) and 
where turbulence is suppressed (thermally dominated mo- 
tion). In this case, any deviations from the expected zero 
mean value of the radial velocity difference between rota- 
tional and inversion molecular transitions in Eq. (|S|) could 
be ascribed to the non-zero A/j, value. 
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In practice, molecular cores are not ideal spheres and 
when observed at higher angular resolutions they frequently 
exhibit complex substructures. The line profiles may be 
asymmetric because of non-thermal bulk motions. Taking 
this into account, the following criteria were formulated: 

1. The line profiles are symmetric described well by a 
single-component Gaussian model (i.e., the minimum 
value of xt ~ 1)- This selection increases the accuracy 
of the line center measurement. Multiple line compo- 
nents may shift the line barycenter and affect the ve- 
locity difference between molecular transitions because, 
e.g., the ratio NH3/HC3N can vary from one component 
to another. 

2. The line widths do not greatly exceed the Doppler width 
because of the thermal motion of material, i.e., the non- 
thermal component (infall, outflow, tidal flow, turbu- 
lence) does not dominate the line broadening. This en- 
sures that selected molecular lines correspond to the 
same kinetic temperature and arise cospatially. For the 
molecules in question we require that the ratio of the 
Doppler 6-parameters, /3 = fe(NH 3 )/&(HC 3 N) or /3 = 
6(NH 3 )/6(N 2 H+), be /3 > 1. 

3. The spectral lines are sufficiently narrow (b ~ 0.1 — 0.2 
km s _1 ) for hyperfine structure components to be re- 
solved. This allows us to validate the measured radial 
velocity by means of different hfs lines of the same 
molecular transition. 

4. The spectral lines are not heavily saturated and their 
profiles are not affected by optical depth effects. The 
total optical depth of the NH3 hf transitions is r < 10. 

We now consider molecular lines that fulfill these ad- 
ditional selection criteria. Among the sources listed in 
Table [TJ two and twelve molecular cores can be selected 
from, respectively, the Medicina and Effelsberg observa- 
tions (Table [5]) , and nine clouds from the Nobeyama data 
set (TableEJ. The sample size is now n — 23, i.e., two times 
smaller than the total data set. 

Table O lists the values obtained from the analysis of 
the NH3 (1,1) and HC3N (2-1) transitions, whereas in 
Table [6] the NH 3 (1,1) and N 2 H+ (1-0) measurements are 
presented. Where individual hfs transitions can be analyzed 
separately, the results of these analyse are also given. The 
data obtained for NH3 are divided formally into 4 groups, 
which are marked in Fig. [JJ and Table by the following 
numbers: 'outer' - 1 and 5, 'inner' -2 and 4, 'central' - 
3, and 'total' combines all hfs lines (Cols. 2-5 of Table [5j 
respectively). For HC3N, the measurements are presented 
in 3 groups, which are marked in Fig. [JJ and Table [3] as: 
'low' - 1, 5, and 6, 'high' - 3 and 4, and 'total' -1, 3, 4, 5, 
and 6 (Cols. 6-8 of Table [SJ respectively). Line number 2 
(F = 1 — 2) is very weak and was never detected in our ob- 
servations. For N2H + , we divided the hfs lines into 3 groups 
on the basis of their relative theoretical intensities. These 
groups are indicated in Fig. [JJ and Tabled as: 'low' - 1, 2, 

5. and 7, 'high' - 3, 4, and 6, and 'total' combines all hfs 
lines (Cols. 6-8 of Table [51 respectively). 

The analyzed molecular line profiles are shown in 
Figs. [2][9l The smooth curves are synthetic spectra calcu- 
lated in the simultaneous fit of all hyperfine components to 
the observed profiles. Bold horizontal lines mark the spec- 
tral ranges included in the ^-minimization procedure. The 
quality of individual fittings can be characterized by the 
normalized xt values given in Tables [5] and [5] for each group 



of measurements. The residuals 'observed data - model' are 
depicted beneath each spectrum in Figs. [2][9] The signal- 
to-noise ratio (S/N) per resolution element shown in these 
figures is calculated at the maximum intensity peak. 

The velocity offsets, AV, and their statistical errors 
are listed in Col. 9 of Tables [5] and [5] When calculat- 
ing AV, we used the mean Vi sr radial velocities based on 
the simultaneous fitting of all hyperfine transitions. The 
AV values estimated in 23 measurements are marked in 
Fig. [JJJ by the circles with la error bars. The filled cir- 
cles indicate sources with thermally dominated motions. 
We also show the ratios of the Doppler 6-parameters: — 
6(NH 3 )/6(HC 3 N) or f3 = 6(NH 3 )/6(N 2 H+), in parentheses 
provide their la errors. For the 32-m and 100-m telescopes, 
AV = V/ sr (HC 3 N) -V/ sr (NH 3 ), whereas for the 45-m tele- 
scope AV = y /sr (N 2 H+) -y /sr (NH 3 ). 

Among the selected clouds, we discover that only one 
(L26O-NH3) violates the selection rule No. 4. The total op- 
tical depth in this case is r = 15.6. We used this cloud since 
the spectrum of NH3 was of high S /N and other selection 
criteria were fulfilled. We note, however, that this cloud has 
the minimum Ay value consistent with a zero offset (see 
Table E] and Fig. [TTJ . 

Both data sets exhibit an excess of positive velocity 
offsets. The most accurate results are obtained for the 
L1498 and L1512 molecular cores observed with the 100-m 
Effelsberg telescope, AV = 21.0 ± 1.5 m s _1 . In this case, 
the systematic error due to the rest frequency uncertainties 
does not exceed 3 m s _1 . The maximum spread between the 
individual hfs Vi sr velocities for NH3 is 5.5 m s _1 (L1498) 
and 4.5 m s _1 (L1512), but it is only 1.0 m s _1 for the 
HC3N hyperfine transitions in both sources. The molecular 
lines are narrow with 6(NH 3 ) = 118 ± 1 m s" 1 , 6(HC 3 N) 
= 91 ± 5 m s" 1 (L1498), and 6(NH 3 ) = 113 ± 1 m s" 1 , 
&(HC 3 N) = 79 ±4 m s^ 1 (L1512), which is in line with the 
assumption that in these two clouds both molecules trace 
the same volume elements. 

The sources with thermally dominated motions (n = 7) 
marked by the filled circles in Fig. [TTJ give the weighted 
mean of (AV) W = 21.1 ± 1.3 m s _1 , the scale 3.4 m s _1 , 
and the median 22 m s -1 . The corresponding M-estimate 
is (AV) M = 21.2 ± 1.8 m s" 1 (scale 4.8 m s" 1 ). 

For the ensemble of n = 23 AV values, we found the 
weighted mean of (AV) W = 20.7±3.0 m s _1 , the scale 14.4 
m s _1 , and the median 22 m s -1 . The robust M-estimate 
is (AV) M = 21.5 ± 2.8 m s" 1 (scale 13.4 m s" 1 ). Thus, in 
the reduced sample, we have good agreement between all 
three estimates of the mean. The scatter in the points (the 
Doppler noise) is lower by a factor of two than that of the 
n = 55 data set. 

The individual data from the 100-m and 45-m telescopes 
provide the following estimates. Effelsberg: (AV) W = 23. 0± 
3.1 m s" 1 (n — 12, scale 10.7 m s~ x , median 22 m s~ x ), 
(AV) M = 23.2 ± 3.8 m s" 1 (scale 13.3 m s" 1 ). Nobeyama: 
(AV) W = 14.0±7.2 m s- 1 (n = 9, scale 21.6 m s^ 1 , median 
22 m s" 1 ), (AV) M = 22.9 ± 4.2 m s" 1 (scale 12.7 m s" 1 ). 

Although the robust M-estimates of the mean for both 
the Effelsberg and Nobeyama observations are consistent, 
the latter has a larger systematic error due to its lower ac- 
curacy at the rest frequencies of the N 2 H + (1-0) transition. 
Taking into account that the rest frequencies of the HC3N 
(2-1) are known with a sufficiently high accuracy (the un- 
certainties in the laboratory and observational frequencies 
are comparable), we take the Effelsberg robust mean as 
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a final value for the velocity offset between the rotational 
and inversion transitions. Being interpreted in terms of the 
electron-to-proton mass ratio variation, this provides the 
value A/i//i= 22 ± 4 stat ± 3 sys ppb (1 ppb = 10~ 9 ). 

5.1. Data reproducibility 

Taking into account numerous perturbation effects dis- 
cussed above and variations in specific parameters such as 
spectral resolution and signal-to-noise ratio, we question 
the consistency of the measured AV values obtained for 
independent telescope systems. This consistency was first 
tested at the 32-m Medicina telescope, where we observed 
two cores, L183 and TMC-1C, in lines of NH 3 (1,1) and 
HC 3 N (2-1) with the ARCOS and MSpecO digital spectrom- 
eters. For L183, we found A7 = 6±7ms" 1 (MSpecO) and 
-4 ±15 m s- 1 (ARCOS), whereas for TMC-1C, the corre- 
sponding quantities were measured to be 67.5 ± 5.8 m s _1 
and 68.5 ± 3.2 m s" 1 . Both results are in good agreement 
to within the la uncertainty interval. 

We also tested the reliability of the velocity offsets by 
obtaining observations of the same cores at different tele- 
scopes. For instance, a quiescent low-mass molecular core 
L1512 in the Taurus Cloud was observed in the NH3 (1,1) 
and HC3N (2-1) lines at the 32-m and 100-m telescopes. 
The corresponding velocity offsets are AV = 17.5 ± 4.1 
m s _1 and 20.5 ± 1.2 m s _1 . Other examples can be found 
in Fig. [TUJ where velocity offsets AV of molecular pairs 
consistent with a singe-component Gaussian model are de- 
picted. 

Some of our targets were partially mapped as indicated 
in Fig.[TU] The six points of TMC-1C are scattered between 
AV = — 7±4 m s~ x and 104±2 m s _1 , which is indicative 
of large systematic shifts caused by bulk motions. Indeed, 
the line widths of NH3 and HC3N demonstrate the dom- 
inant influence of the non-thermal component. The mea- 
sured ratio (3 — 6(NH3)/6(HC3N) at these points is equal to 
1.0, 0.8, 1.2, 1.0, 1.2, and 1.0 at, respectively, (Aa, AS) = 
(0,0), (0,-40), (0,40), (-120, 90), (-60, 45), and (60,-45) 
arcsec. Similar velocity offsets are observed for L1495 and 
L183. In general, the scatter for the whole sample covers the 
range —100 < AV^ < 100 m s _1 . However, the scatter for 
the subsample selected in accord with the additional selec- 
tion criteria from Sect. [5] decreases considerably. The corre- 
sponding points shown in Fig. [TT] are distributed across the 
interval from — 10 to 43 m s _1 . Thus, as expected, molec- 
ular cores with thermally dominated gas motions are the 
most suitable targets for precise measurements of Afi/ p, . 

It is essential that the observed spectral-line displace- 
ments are based on differential measurements of radial ve- 
locities recorded at the same telescope over short periods 
of time. It is also important that different receivers used in 
these observations enable the independent Doppler tracking 
of the observed molecular lines. Thus, there is no need for 
an accurate definition of the zero point, and only internal 
precisions of the measured line positions and the uncertain- 
ties in the rest- frame frequencies restrict the accuracy of the 
measured velocity offset. 

5.2. Comparison with previously obtained results 

In our preliminary analysis (LMK), we used high quality ra- 
dio spectra of molecular cores in lines of NH3 (J, K) = (1,1) 



(23 GHz), CCS J N = 2i - 1 (22 GHz), HC 3 N J = 5-4 (45 
GHz), and N 2 H+ J = 1-0 (93 GHz). The observations were 
carried out with the 100-m Green Bank Telescope (GBT) 
by Rosolowsky et al. (2008) for the Perseus molecular cloud, 
by Rathborne et al. (2008) for the Pipe Nebula, and with 
the Nobeyama 45-m telescope by Sakai et al. (2008) for 
infrared dark clouds (IRDCs). 

The most accurate estimates were obtained from care- 
fully selected subsamples of the NH 3 /CCS pairs observed 
with high spectral resolution (FWHM = 25 m s" 1 ) in the 
Perseus molecular cloud (n = 21) and the Pipe Nebula 
(n = 8): AV = 52±7 st at±14 S y S m s" 1 and 69±ll stat ±14 sys 
m s _1 , respectively. The analysis of the n — 36 NH3/N2H+ 
pairs and n — 27 NH3/HC3N pairs observed at lower 
spectral resolution (FWHM = 120-500 m s" 1 ) in the 
IRDCs inferred that AV = 160 ± 32 sta t ± 14 sys m s" 1 
and 120 ± 37 sta t ± 31 sys m s _1 , respectively. At the same 
time, the rotational-rotational velocity differences for the 
IRDCs did not show any significant offset: AV = F(N2H + ) 

- ^(HC 3 N) = -17 ± 34 km s" 1 . 

The physical properties of low- mass (M < 10M Q ) 
molecular cores observed with the GBT are similar to those 
analyzed in the present paper. However, the IRDCs are 
more massive objects (M > 100Af Q ) with high turbulent 
velocities (^1-3 km s _1 ). The IRDCs are also more dis- 
tant clouds, which are located towards the Galactic cen- 
ter at distances of 2-5 kpc. The spectral resolution used 
in observations of the IRDCs does not allow us to resolve 
the hyperfine structure of the HC3N J = 5 — 4 transition 
(Sakai et al. 2008). The spectra of the HC 3 N lines were 
fitted to a single Gaussian function and the corresponding 
radial velocities were calculated based on a single hyperfine 
line J = 5-4,F = 5- 4 with the frequency 45490.3102(3) 
MHz taken from the JPL catalog^. As a result, a systematic 
shift related to unresolved hyperfine components was intro- 
duced. We also note that the JPL frequencies are systemati- 
cally shifted with respect to the CDMS catalogue (Miiller et 
al. 2005). The CDMS frequency of the J = 5 - 4, F = 5 - 4 
transition is, for example, 45490.3137(5) MHz. Thus, it 
seems plausible that the large offset, AV > 100 m s _1 , 
deduced from the IRDCs is caused by the superposition of 
unaccounted systematic shifts. 

As for the low-mass cores in the Perseus molecular cloud 
and the Pipe Nebula, the discrepancy between our present 
estimate, AV^ ~ 20 m s _1 , and the previously obtained 
value, AV ~ 60 m s" 1 , is probably caused by uncertain- 
ties in the rest-frame frequency of the CCS Jm = 1\ — Iq 
transition. In particular, the frequency 22344.033(1) MHz 
of this transition (used by both Rosolowsky et al. and 
Rathborne et al.) was calculated by comparing with the 
HC 3 N J = 5 — 4 line observed towards a cold dark cloud 
L1498 assuming that (i) CCS and HC3N are co-spatially 
distributed and (ii) the rest-frame frequency of HC3N 
J = 5 - 4, F = 5 - 4 is 45490.316(1) MHz (Yamamoto 
et al. 1990). If the HC3N frequency is 2 kHz lower as in- 
dicated in the CDMS catalogue, then the laboratory fre- 
quency of CCS decreases by 1 kHz, where 1 kHz at 22.3 
GHz corresponds to 13.4 m s _1 , i.e., the difference F(CCS) 

- F(NH3) becomes equal to ~ 46 m s _1 . The only known 
laboratory measurement of the CCS frequency infers, how- 
ever, that 22344.029(4) MHz (Lovas et al. 1992). With this 
value V^CCS) - V^NHs) ~ 6 m s^ 1 , but its error reaches 

6 http://spec.jpl.nasa.gov 
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53 m s 1 (4 kHz at 22.3 GHz). The rest-frame frequency of 
CCS reported in the CDMS is 22344.0308(10) MHz. With 
this frequency, we obtain for the Perseus molecular cloud 
the offset AV = 22 ± 7 stat ± 14 sys m s" 1 , which is in full 
agreement with our present estimate, AV = 23±4 sta t±3 sys 
m s _1 . This consideration makes it quite obvious that to in- 
crease the reliability of the results obtained, high precision 
laboratory measurements of molecular transition frequen- 
cies are badly needed. 

6. Conclusions and future prospects 

Using the 32-m Mcdicina, 100-m Effclsberg, and 45-m 
Nobeyama radio telescopes, we have performed precise 
spectral measurements of the relative radial velocity dif- 
ferences between the rotational transitions in HC3N (J = 
2 — 1) and N2H 4 " (J = 1—0) and the inversion transi- 
tion in NH3 (J,K) = (1, 1). We detect a velocity offset of 
AV — 23 ± 4 sta t ± 3 sys m s _1 between the rotational and 
inversion lines observed in cold molecular cores with dom- 
inating thermal motion. We do not find any plausible sys- 
tematic effects that could mimic an offset of about 20 m s _1 
between rotational and inversion transitions and would be 
regularly reproduced in observations of different cold molec- 
ular cores with different facilities. The measured positive 
offset is qualitatively consistent with our preliminary re- 
sult, AV ~ 60 m s" 1 , obtained on the basis of the GBT 
observations of Rosolowsky et al. (2008) and Rathborne 
et al. (2008) of the NH3/CCS pairs in molecular cores in 
the Perseus molecular cloud and the Pipe Nebula (LMK). 
However, the rest-frame frequency of the CCS Jn = %i — lo 
transition is not well known, and before deciding whether 
the difference between the current and preliminary results 
is statistically significant or not, new high precision labora- 
tory measurements of the CCS frequency should be carried 
out. 

If we assume that the measured velocity offset is 
caused by the electron-to-proton mass ratio variation, then 
A^i I ii~ 22 ± 4 sta t ± 3 sys ppb. To account for the conditions 
of the terrestrial laboratory experiments, the non-zero A/i 
would require chameleon-like scalar field models that pre- 
dict a strong dependence of mass and coupling constant on 
the ambient matter density. 

In any case, our estimate defines a strongest conser- 
vative upper limit to |A/i//i| < 3 x 10~ 8 , which can be 
considered as a reference point at z = 0. We note that ex- 
tragalactic molecular clouds have gas densities similar to 
those in the interstellar clouds of the Milky Way. Thus, 
the value of A{i//i in high-z molecular systems is expected 
to be at the same level as in the interstellar clouds, i.e., 
A/i//i^ 10~ 8 , provided that no temporal dependence of 
the electron-to-proton mass ratio is present. 

To be completely confident that the derived velocity 
shift is not caused by kinematic effects in the clouds but 
reflects the density-modulated variation of A/j,/ fi , new high 
precision radio-astronomical observations are needed for a 
wider range of objects. These observations should also tar- 
get different molecules (i.e., not NH3) with tunneling tran- 
sitions sensitive to the changes in fi and molecules with 
A-doublct lines that also exhibit enhanced sensitivity to 
variations in [i and a (Kozlov 2009). 

For the molecular transitions in question, it is very im- 
portant to measure their rest-frame frequencies with an ac- 
curacy of about 1 m s _1 in laboratory experiments. In some 



Table 2. Hypcrfine components of the NH3 (J, K) = 
(1,1) transition. The velocity offsets, Vi, are relative to 
23.694495487 GHz. Data from Kukolich (1967) and Ho 
(1977). Shown in parentheses are la errors. 



# F[ F' F x F Weight v iy m s" 1 



1 

2 



3 



4 



5 



^The error of 0.35 m s _1 corresponds to the difference 
between theoretically predicted (Ho 1977) and astrono- 
mically determined (Rydbeck et al. 1977) frequencies. 

Table 3. Hyperfine components of the HC3N J = 
2 — 1 transition. The velocity offsets, Vi 1 are relative 
to 18.19621694 GHz. Data from Lapinov (2008, private 
comm.) and Miiller et al. (2005). 



# 


F' 


F 


Weight 


Vi, ms 


-1 


1 


1 


1 


0.333 


-35548.771 


'2.82 


2 


1 


2 


0.022 


-14167.971 


"2.83 


3 


3 


2 


1.867 


-1540.961 


'2.82 


4 


2 


1 


1.000 


O.OOi 


'2.82 


5 


1 





0.444 


17806.601 


'2.82' 


6 


2 


2 


0.333 


21380.641 


"2.82 



cases the present uncertainties in the rest-frame frequencies 
are larger than the errors of radio-astronomical measure- 
ments, thus preventing unambiguous conclusions. 

In addition, searches for variations in the fine-structure 
constant a within the Milky Way disk using mid- and 
far-infrared fine-structure transitions in atoms and ions 
(Kozlov et al. 2008), or searches for variations in the com- 
bination of a 2 / \i using the [C 11] A158 /xm line and CO rota- 
tional lines (Levshakov et al. 2008b), or the [Cl] A609 /xm 
line and low-lying rotational lines of 13 CO (Levshakov et 
al. 2009) would be of great importance for cross-checking 
these results. 
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\ 


1 ^ 


n 
u 


\J.O 


n ^1 7 




1 


0.5 





0.5 


0.159 


-19411.73(61) 


1 


1.5 


2 


2.5 


0.357 


-7815.25(63) 


1 


1.5 


2 


1.5 


0.040 


-7372.80(75) 


1 


0.5 


2 


1.5 


0.198 


-7233.48(58) 


2 


1.5 


2 


2.5 


0.071 


-250.92(60) 


1 


1.5 


1 


1.5 


0.198 


-213.00(47 


2 


2.5 


2 


2.5 


1.000 


-132.38(47 


1 


0.5 


1 


1.5 


0.040 


-75.17(75) 


2 


1.5 


2 


1.5 


0.643 


192.27(56 


2 


2.5 


2 


1.5 


0.071 


311.03(80 


1 


1.5 


1 


0.5 


0.040 


322.04(35 


1 


0.5 


1 


0.5 


0.079 


460.41(46) 


2 


1.5 


1 


1.5 


0.040 


7351.32(49) 


2 


2.5 


1 


1.5 


0.357 


7469.67(73) 


2 


1.5 


1 


0.5 


0.198 


7886.69(72) 





0.5 


1 


1.5 


0.317 


19315.90(71) 





0.5 


1 


0.5 


0.159 


19851.27(35) 
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Table 1. Target list. 



# 


Source 


R.A. 
(J2000.0) 


Dec. 
(J2000.0) 


Vlsr 

(km s- 1 ) 


Obs. 


a 


Molecules 6 


References c 


1 


L1355 


02:53:12.1 


+68:55:52 


-3.7 


m, e 




NH 3 




14 


2 


NH3SRC8 


03:25:26.3 


+30:45:05 


+4.1 


m, e, 


n 


NH 3 , 


N 2 H + 


5 


3 


NH3SRC28 


03:27:26.4 


+29:51:08 


+5.6 


e 




NH 3 




5 


4 


NH3SRC36 


03:27:55.9 


+30:06:18 


+4.7 


m 




NH 3 




5 


5 


JNH3SRC91 


no o r\ 1 r ri 

03:30:15.2 


i on no on 

+30:23:39 


+5.9 


m 




INrl3 




5 


6 


NH3SRC95 


03:30:32.1 


+30:26:19 


+6.1 


m 




INrl3, 




5 


-7 

7 


JNH3SKC97 


no on cn c 

03:30:50.5 


i on ACl IT 

+30:49:17 


+7.7 


e 








5 


8 


L1489 


n a n a a r\ r\ 

04:04:49.0 


i o/" 1 -i o An 

+26:18:42 


+6.7 


e, n 




INrl3, 


TT/' t 1S.T AT T_T — 

rlL,3lN , 1N2J4 


1, 4 


9 


L1498 


04:10:51.4 


+25:09:58 


+7.8 


e, n 




INrl3, 


£103 IN , IN 2-tl 


1, 3, 4, 6 


10 


L1495 


04:14:08.2 


+28:08:16 


+6.8 


e, n 




MTT 

JNH 3 , 


TJf~* "NT M U + 

rlL-3iN , IN 2Ti 


1, 6 


11 


L1521F 


04:28:39.8 


+26:51:35 


+6.5 


e 






rl03lN 


4, 6, 10 


12 


LI 524 


04:29:22.7 


+24:34:58 


+6.5 


e, n 




1VTTJ 


T 1 AT AT TJ — 

rlOsiN , IN21I 


4, 8, 10 


13 


B217S- W 


U4:z7:4o.4 


i 1^7 CO 

+26:17:52 


+7.0 


e 




i\tu 
1NH 3 , 


TTTI AT 

rl03iN 


1, 4, 12 


1 A 

t4 


L14U0K 


n /i .on. co n 


+o4:ol:oo 


+3.3 


e, n 




lNrl 3 , 


TT/' AT 

rl03lN 


1 A a 

1, 3, 4, 6 


10 


1 Mt-z 


U4:oz:4o.o 


+24:2o:12 


1 t! O 

+0.2 


e, n 




lNrl 3 , 


IT/ 1 AT AT TT+ 

rl03iN , IN 2Jl 


t 


16 


L1536 


04:33:23.5 


+22:42:46 


+5.6 


e, n 




lNrl 3 , 


\~X ( ^ AT AT TJ + 

rl03lN , IN 2T1 


1, 4 


17 


CB22 


C\ A A n O O 1 

04:40:33.1 


i on r* r* r\ a 

+29:55:04 


+6.0 


e 




iMn 3 , 


ri v,/ 3 In 


4 


18 


1JV1C-1C 


04:41:38.8 


i r»/° nn a n 

+26:00:42 


+5.2 


m, e, 


n 


i\tu 
1NH 3 , 


T 1 AT AT TJ + 

£103 IN 5 J-N 2-tl 


1, 2, 3, 4 


ty 


1 MIj- t / rlG31N 


C\A . A 1 . A 1 O 

U4:4t:4t.o 


i oc A i . /in 
+2o:4t:42 


i c n 
+0.9 


m 




1N11 3 , 


TTTI AT 

£103 IN 


n 11 
9, 11 


on 
2U 


rinno 

OB2o 




i on. on. 1 1 
+29:39:11 


i tt n 
+0.0 


e, n 




NH 3 . 


TTin AT AT TT - ! - 

HC 3 JN, N2H T 


4, D 


21 


L1517B 


04:55:18.8 


+30:38:04 


+5.8 


e, n 




NH 3 , 


rr / < -iw -r -\t tt_|_ 

HC 3 N, N2H + 


2, 4, 6 


22 


L1512 


05:04:09.6 


+32:43:09 


+7.1 


m, e, 


n 


NH 3 , 


HC 3 N, N 2 H+ 


1, 2, 4, 6, 8 


23 


L1544 


05:04:16.5 


+25:10:48 


+7.2 


e, n 




NH 3 , 


HC 3 N 


3, 4, 6, 8, 9 


24 


LI 34 A 


15:53:36.2 


-04:35:26 


+2.7 


e 




NH 3 




4, 6, 8 


25 


L183 


15:54:06.4 


-02:52:23 


+2.4 


m, e, 


n 


NH 3 , 


HC 3 N, N 2 H+ 


4, 6, 7, 8 


26 


L43 


16:34:35.0 


-15:46:36 


+0.7 


e, n 




NH 3 , 


HC 3 N, N 2 H+ 


4 


27 


L260-NH3 


16:47:06.6 


-09:35:21 


+3.5 


e, n 




NH 3 , 


HC 3 N, N 2 H+ 


1, 8, 12 


28 


L234A 


16:48:06.8 


-10:51:48 


+2.9 


e 




NH 3 , 


HC 3 N 


4,8 


29 


L234S-E 


16:48:08.6 


-10:57:25 


+3.2 


e 




NH 3 




4 


30 


L63 


16:50:15.4 


-18:06:06 


+5.7 


e 




NH 3 , 


HC 3 N 


4, 8 


31 


CB68 


16:57:19.4 


-16:09:25 


+5.2 


e 




NH 3 , 


HC 3 N 


4 


32 


L492 


18:15:49.6 


-03:46:14 


+7.7 


e 




NH 3 , 


HC 3 N 


6 


33 


T 770 

L778 


i n or* oo 7. 

19:2o:32.5 


i no cro /in 

+23:58:42 


+9.9 


m 




NH 3 




1, 4 


34 


B335 


19:37:01.3 


+07:34:07 


+8.4 


e 




NH 3 , 


HC 3 N 


4, 8 


35 


L694-2 


19:41:04.4 


+10:57:02 


+9.5 


m, e 




NH 3 




6, 13 


36 


L1049-2 


20:41:30.5 


+57:33:44 


+0.1 


e 




NH 3 




15 


37 


L1251T3 


22:29:50.3 


+75:13:25 


-3.8 


e, n 




NH 3 , 


HC 3 N, N 2 H+ 


4 


38 


L1251T2 


22:30:31.2 


+75:13:38 


-4.0 


e, n 




NH 3 , 


HC 3 N, N 2 H+ 


4 


39 


L1251T1 


22:31:02.3 


+75:13:39 


-4.2 


e, n 




NH 3 , 


HC3N, N 2 H+ 


4 


40 


L1251C 


22:35:53.6 


+75:18:55 


-4.7 


e, n 




NH 3 , 


HC 3 N, N 2 H+ 


16,17 


41 


L1197 


22:37:02.1 


+58:57:22 


-3.1 


e, n 




NH 3 , 


HC3N, N 2 H+ 


6,13 



a Observations: m - Medicina 32-m, e - Effelsberg 100-m, and n - Nobeyama 45-m radio telescopes. 
b Detected transitions: NH 3 (J, K) = (1, 1), HC3N J = 2 - 1, N 2 H+ J = 1 — 0. 

c References: (1) Benson & Myers 1989; (2) Myers et al. 1983; (3) Howe et al. 1994; (4) Jijina et al. 1999; 
(5) Rosolowsky et al. 2008; (6) Crapsi et al. 2005; (7) Pagani et al. 2009; (8) Fuller & Myers 1993; 
(9) Churchwell et al. 1984; (10) Codella et al. 1997; (11) Benson & Myers 1983; (12) Motte & Andre 2001; 
(13) Lee et al. 2007; (14) Park et al. 2004; (15) Lee et al. 1999; (16) Caselli et al. 2002; (17) Walsh et al. 2004. 

Appendix A: Robust statistics 

A comprehensive review of robust statistical procedures is 
given by Hampel et al. (1986) (see also Chap. 15.7 in Press 
et al. 1992). 

Robust estimates are insensitive to departures of the 
distribution shape from the assumed data distribution. In 
particular, M-estimates (maximum-likelihood) are usually 
relevant for the estimate of parameters with continuous dis- 
tributions. The shift (mean), m, and the scale (standard de- 
viation), s, of a sample of data points {xi} are estimated by 
minimizing the sum ^[(x^— m)/s]. Here ^ is a weighting 
function such that the weights given to individual points 



first increase with deviations from m, and then decrease 
so that very outlying points (outliers) are not counted at 
all (re-descending M-estimate). The following functions are 
prescribed for use as >]/: 

Hampcl's function 



*(u) = sign(u) { 8 ; 5 1| M | 



u\, 0<M<1.7 

1.7, 1-7 < |uj < 3.4 

. 3.4 < M < 8.5 

0. " u > 8.5 
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NH, (1,1) 



Main Lines, At'-O 



Satellite Lines, 



Sate] :te lines 



'S -20 -15 -10 -5 



10 15 



I 

1 HC 3 N (2-1) 


1 i i i i | i i i i | 


3-2 














8-1 




7F=1-1 










1° T 




1 2 










-1 

r , , 


8 


3 4 


5 6 : 



-36 -30 -25 -20 -15 -10 -5 5 10 15 20 



N,H* (1-0) 



Velocity Shift (km s" 1 ) 



Fig. 1. Sketch of the hyperfine transitions for the NH3 (1,1), HC3N (2 — 1), and N2H+ (1 — 0) states. Line identification 
numbers are given in each panel (see Tables [2] - [4]) . 



Table 4. Hyperfine components of the N2H + J = 
1 — transition. The velocity offsets, Vi, are relative to 
93.1737722 GHz. Data from Miiller et al. (2005). 



# 


F[ 


F' 


Fi 


F 


Weight 


Vi, m s 1 


1 





1 


1 





0.076 


-8008.2(13.8) 







1 


1 


1 


0.100 


-8008.2(13.8) 







1 


1 


2 


0.252 


-8008.2(13.8) 


2 


2 


1 


1 


2 


0.017 


-622.6(13.8) 




2 


1 


1 





0.134 


-622.6(13.8) 




2 


1 


1 


1 


0.278 


-622.6(13.8) 


3 


2 


3 


1 


2 


1.000 


0.00(13.2) 


4 


2 


2 


1 


2 


0.111 


956.9(13.5) 




2 


2 


1 


1 


0.604 


956.9(13.5) 


5 


1 


1 


1 


1 


0.051 


5541.6(13.5) 




1 


1 


1 


2 


0.159 


5541.6(13.5) 




1 


1 


1 





0.218 


5541.6(13.5) 


6 


1 


2 


1 


1 


0.111 


5982.4(13.5) 




1 


2 


1 


2 


0.604 


5982.4(13.5) 


7 


1 





1 


1 


0.143 


6934.2(13.2) 



Andrew's sine 

*(u) = /™(c)' H <CT 

W 10, M > C7T 

optimal value for c = 2.1, and 
Tukey's biweight 

*H = { U ( 1 -^) 2 ' H <c 

[0, \u\>c 



optimal value for c = 5 — 6. 

The minimization occurs iteratively, the initial guess 
for to being the median of the data set, and for the scale 
s the normalized median absolute deviation (1.48-MAD). 
All estimations presented in the paper were obtained with 
Tukey's biweight. Other ^-functions provide similar results. 
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Fig. 7. Spectra of NH 3 (1,1) and N 2 H+ (1 - 0) toward the cores L1498, L1536, and CB23 obtained at the Nobeyama 
45-m radio telescope. The histogram shows the data, the solid curve shows the fit, and the residual is plotted below 
each profile. The horizontal thick bars mark spectral windows used in the fitting procedure. The data are the arithmetic 
means of all the observations. The size of the resolution element (pixel) is 49 m s _1 for NH3, and 25 m s _1 for N2H" 1 ". 
For each spectrum, the signal-to-noise ratio (S/N) per pixel at the maximum intensity peak is depicted. 
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Table 5. Radial velocities, Vi sr , Doppler parameters, b, and corresponding xt values normalized per degree of freedom. 
Data from the Medicina 32-m and Effelsberg 100-m radio telescopes. The numbers in parentheses correspond to ler 
statistical errors. 
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Table 6. Radial velocities, Vi sr , Doppler parameters, 6, and corresponding xl values normalized per degree of freedom. 
Data from the Nobeyama 45-m radio telescope. The numbers in parentheses correspond to la statistical errors. 
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0.092(15) 

0.84 

5.629(11 
0.121(24 
0.76 
6.0515(46) 
0.097(13) 
0.64 

5.8190(50) 

0.109(21) 

1.00 

7.1290(45 
0.0847(92 
0.89 

2.4395(50) 

0.149(17) 

0.91 

2.4190(26 
0.1325(64 
1.18 

0.7485(31 
0.1368(88 
1.81 

3.4700(59 
0.0997(82 
1.34 



7.8400(43 
0.0987 40 
1.16 

5.6295(92) 

0.117(40) 

0.77 

6.0385(45) 

0.119(12) 

1.40 

5.8215(53) 

0.138(16) 

1.03 

7.1260(41) 

0.085(11) 

0.89 

2.4485(45 
0.1484 61 
0.82 

2.4135(26) 
0.1624(69) 
1.01 

0.7470(25) 

0.173(17) 

0.90 

3.4710(52) 

0.111(26) 

1.12 



7.8430(30 
0.1025(68 
1.08 

5.6290(80) 

0.118(20) 

0.75 

6.0460(33 
0.1043(92 
1.00 

5.8195(38) 

0.127(14) 

1.03 

7.1270(35) 
0.0850(90) 
0.86 

2.4445(32) 

0.154(56) 

0.91 

2.4165(18 
0.1423 54 
1.30 

0.7480(19) 
0.1512(93) 
1.60 

3.4690(37) 
0.1061 79) 
1.24 



+23.0(3.4) 
+33.5(8.1) 
+11.5(3.6) 
+21.0(5.1) 
+24.5(4.3) 
+22.0(3.6) 
-10.0(2.1) 
+43.0(2.5) 
-10.0(4.0) 



a Offset (Aa, A<5) = (-31.5", 25"). 



